Electronic Version 
Stylesheet Version vl.1.1 



Description 



[A METHOD FOR SUPPLY CHAIN 
DECOMPOSITION] 

Cross Reference to Related Applications 

[0001] The present application is related to pending U.S. Patent 
Application 10/ , filed concurrently herewith to Den- 
ton et al., entitled " A Method for Supply Chain Compres- 
sion" having (IBM) Docket No. BUR920030197US1; U.S. 

Patent Application 10/ , filed concurrently herewith 

to Denton et al., entitled "METHOD FOR PURCHASE ORDER 
RESCHEDULING IN A LINEAR PROGRAM" having (IBM) 
Docket No. BUR92004009US1; U.S. Patent Application 

10/ , filed concurrently herewith to Denton et al., 

entitled "A METHOD FOR OPTIMIZING FOUNDRY CAPAC- 
ITY" having (IBM) Docket No. BUR920030195US1; U.S. 

Patent Application 10/ , filed concurrently herewith 

to Denton et al., entitled "METHOD FOR FAIR SHARING 
LIMITED RESOURCES BETWEEN MULTIPLE CUSTOMERS" 
having (IBM) Docket No. BUR920040010US1; U.S. Patent 



Application 10/ , filed concurrently herewith to Den- 
ton et al., entitled "A METHOD FOR CONSIDERING HIERAR- 
CHICAL PREEMPTIVE DEMAND PRIORITIES IN A SUPPLY 
CHAIN OPTIMIZATION MODEL" having (IBM) Docket No. 

BUR920030198US1; U.S. Patent Application 10/ , 

filed concurrently herewith to Denton et al., entitled " 
Method for Simultaneously Considering Customer Commit 
Dates and Customer Request Dates" having (IBM) Docket 
No. BUR920040008US1; and U.S. Patent Application 

10/ , filed concurrently herewith to Denton et al., 

entitled "METHOD FOR IDENTIFYING PRODUCT ASSETS IN 
A SUPPLY CHAIN USED TO SATISFY MULTIPLE CUSTOMER 
DEMANDS" having Docket No. BUR820030346US1. The 
foregoing applications are assigned to the present as- 
signee, and are all incorporated herein by reference. 
Background of Invention 

[0002] p ie id 0 f tne invention 

[0003] The present invention relates to computer implementable 
decision support systems for determining an optimal sup- 
ply chain plan. General methodologies within this field of 
study include advanced planning systems, optimization 
and heuristic based algorithms, constraint based pro- 



gramming, and simulation. 

[0004] Description of the Related Art 

[0005] a fundamental problem faced in all manufacturing indus- 
tries is the allocation of material and capacity assets to 
meet end customer demand. Production lead times neces- 
sitate the advance planning of production starts, inter- 
plant shipments, and material substitutions throughout 
the supply chain so that these decisions are coordinated 
with the end customers" demand for any of a wide range 
of finished products (typically on the order of thousands 
in semiconductor manufacturing). Such advance planning 
depends upon the availability of finite resources which in- 
clude: finished goods inventory, work in process inventory 
(WIP) at various stages of the manufacturing system, and 
work-center capacity. Often, there are alternative possi- 
bilities for satisfying the demand. Products may be built at 
alternative locations and within a location there may be 
choices as to which materials or capacity to use to build 
the product. The product may be built directly or acquired 
through material substitution or purchase. When limited 
resources prevent the satisfaction of all demands, deci- 
sions need to be made as to which demand to satisfy and 
how to satisfy it. This resource allocation problem is often 



addressed through linear programming. 
[0006] The below-referenced U.S. Patents disclose embodiments 
that were satisfactory for the purposes for which they 
were intended. The disclosures of both the below- 
referenced prior U.S. Patents, in their entireties, are 
hereby expressly incorporated by reference into the 
present invention for purposes including, but not limited 
to, indicating the background of the present invention and 
illustrating the state of the art: U.S. Patent 5,971,585, 
"Best can do matching of assets with demand in micro- 
electronics manufacturing," October 26, 1999; U.S. Patent 
5,943,484, "Advanced material requirements planning in 
microelectronics manufacturing," August 24, 1999; and 
Nemhauser, G.L. and Wolsey, LA., 1999, Wiley, Integer and 
Combinatorial Optimization . 
Summary of Invention 

[0007] The invention provides a method and system for solving a 
linear program having constraints in a production plan- 
ning system. The invention first determines which of the 
constraints can be temporarily relaxed based on stocking 
point criteria. The stocking point criteria relates to time 
dependent stocking points that include part numbers, lo- 
cations of parts identified by the part numbers, and the 



time periods when the parts will be available. The inven- 
tion relaxes the constraints that can be relaxed and de- 
composes the linear program into smaller independent 
linear programs. The invention initially solves the smaller 
independent linear programs with relaxed constraints 
(simultaneously in parallel) to produce an initial solution. 
Next, the invention replaces variables in the linear pro- 
gram with constants based on this initial solution. After 
this, the invention restores the material balance and 
sourcing constraints and finally solves (re-solves) the lin- 
ear program using the constants and with all constraints 
in place to obtain a complete solution of the linear pro- 
gram. 

[0008] jhe process of determining which of the constraints can 
be temporarily relaxed identifies constraints that are as- 
sociated with the least complex parts within bills— 
of-materials used by the linear program as the constraints 
that can be temporarily relaxed. Additionally, this process 
identifies constraints that are associated with parts that 
have supply availability and lack capacity constraints as 
the constraints that can be temporarily relaxed. Further, 
this process identifies constraints that are associated with 
parts that are available during the planning horizon of the 



linear program as the constraints that can be temporarily 
relaxed (the planning horizon includes an initial planning 
horizon, shipping lead time, and manufacturing cycle 
time). 

Brief Description of Drawings 

[0009] Figure 1: Overview of the structure of a linear program- 
ming application. 

[0010] Figure 2: Illustration of relaxation of material balance 
equations at the bottom of the bill-of-material and re- 
sulting potential for decomposition. 

[001 1] Figure 3: Illustration of the flow of the present invention. 

Detailed Description 

[0012] Detailed Description of Preferred Embodiment(s) of the In- 
ventionThe invention solves a supply-chain optimization 
model according to a two stage algorithm. Initially, a sub- 
set of material balance and sourcing constraint equations 
are relaxed to reduce the size of the linear program and 
also to increase the likelihood that it can be decomposed 
into multiple smaller linear programs. Next, the relaxed 
linear program is solved (via decomposition if possible) 
using an existing methodology (e.g., simplex method). 
The current solution is fixed, the previously relaxed con- 



straints (and associated variables) are added back to the 
model, and the remaining (much smaller) model is solved 
to determine the final solution. 

[0013] jo contrast the present invention, a conventional produc- 
tion planning linear program "LP" is shown below (such as 
that described in U.S. Patent 5,971,585, which is incorpo- 
rated herein by reference). This LP makes decisions in- 
cluding: production starts, material substitutions, and 
shipments planned to customers, between manufacturing 
and distribution locations, and from vendor suppliers. A 
LP is composed of an objective function that defines a 
measure of the quality of a given solution, and a set of 
linear constraints. The types of equations used in produc- 
tion planning models are well know to those practiced in 
the art and include:(l) Material Balance Constraints, which 
ensure conservation of material flow through the network 
of stocking points comprising the supply chain. 

[0014] (2) Capacity Constraints, which ensure that the capacity 
available for manufacturing activities is not exceeded. 

[0015] (3) Backorder Conservation Constraints, which balance the 
quantity of a given part backordered in a given planning 
period with the quantity backordered in the previous 
planning period and the net of new demand and new 



shipments. 

[0016] (4) Sourcing Constraints, which define target ranges 

(minimum and maximum) of shipments that should be 
made from a particular manufacturing or vendor location 
in the supply chain. 

[0017] a conventional LP formulation is provided below in the 
form familiar to those practiced in the art; i.e., definition 
of subscripts, definition of objective function coefficients, 
definition of constants, definition of decision variables, LP 
formulation or equations. 

[001 8] Definition of Subscripts 

[0019] j_ time period 

[0020] m _ material (part number) 

[0021] a _ p | ant location within the enterprise 

[0022] n material being substituted 

[0023] z _ group (which represents a family or collection of part 
numbers) 

[0024] e process (a method of purchasing or manufacturing a 
material at a plant)v- receiving plant location 

[0025] £ demand center (i.e., customer location) (Note: the set of 
customer locations is mutuallyexclusive from the set of 



plant locations) 
[0026] q _ demand class which indicates relative priority 

[0027] w _ resource capacity which could be a machine, labor 
hour, or other constraint 

[0028] u represents a consumer location which refers to an inter- 
nal plant, external demand center, or toa generic indicator 
meaning any plant/or demand center 

[0029] Definition of Objective Function Coefficients 

[0030] prc. - cost of releasing one piece of part m during pe- 

jmae 

riod jat plant a using process e 
[0031] SUBC substitution cost per piece of part number n 

jmna - 

which is being substituted by part number m during pe- 
riodjat plant a 

[0032] jq _ transportation cost per piece of part number m 

jmav 

leaving plant a during period j which are destined for 
plant v 

[0033] invc - inventory cost of holding one piece of part 

jma * 

number m at the end of period j at a particular plant a 
[0034] DMAXC - cost per piece of exceeding the maximum 

jzau 

amount of shipments of group z parts from plant a to 
consuming location(s) u during period j 
[0035] DMINC - cost per piece of falling short of the minimum 

jzau 



amount of shipments specified for group z parts from 
plant a to consuming location(s) u during period j 
[0036] BOC - backorder cost of one piece of part m at the end 

jmkq 

of period j for class q demand at customer location k 

[0037] Definition of Constants 

[0038] DEMAND - demand requested during time period jfor 

jmkq 

part number m at customer location k for demand class a 
[0039] RECEIPT - quantity of projected wip and purchase order 

jma 

receipts for part number m expected to be received at 
plant a during time periodj 

[0040] CAPACITY - Capacity of resource w available at plant a 

jaw 

during period jto support production starts 
[0041] CAPREQ. - Capacity of resource w required for part 

number m at plant a for process e during period j 
[0042] QTYPER - quantity of component m needed per part 

jmaen 

number n during period j at plant a using process e 
[0043] YIELD - output of part number m per piece released or 

jmae 

started at plant a during time period j using process e 
[0044] SUBQTY - quantity of part number m required to sub- 

jmna 

stitute for one piece of part number n at plant a during 
time period j 

[0045] MAXPCT - maximum percentage of total shipments of 

jzau 

group z (collection of parts) leaving supplier a during pe- 



riod jto support consumption at consuming location(s) u 
[0046] MINPCT - minimum percentage of total shipments of 

jzau 

group z (collection of parts) leaving supplier a during pe- 
riod jto support consumption at consuming location(s) u 
[0047] qj _ Cycle time. The number of periods between the 

jmae 

release and completion of part m jobs for releases made 
using process e at plant a during time period j 
[0048] jj _ transport time for part number m from plant a to 

mav 

plant v 

[0049] Definition of LP Decision Variables 

[0050] | _ inventory at the end of period j for part number m at 

jma 

a particular plant a 
[0051] p _ Production starts of part m during period jat plant 

jmae 

a using process e 
[0052] |_ _ Quantity of part number n which is being substi- 

jmna 

tuted by part number m during periodj at plant a 
[0053] j _ internal shipments of part number m leaving plant 

jmav 

a during period j which are destined for plant v 
[0054] p - Shipments of part number m leaving plant a dur- 

jmakq 

ing period jand satisfying class g demand at external cus- 
tomer k 

[0055] b _ Back orders of part m at the end of period j for 

jmkq 

class a demand at customer location k 



[0056] h _ Total shipments of group z (z is a "collectiorfof 

jzu 

parts) leaving suppliers during period j to support con- 
sumption at consuming location(s) u 
[0057] s - Amount by which total shipments of parts in z from 

jzau 

plant a to consuming location(s) u during period j exceeds 
the maximum amount specified as desired in the sourcing 
rules 

[0058] q _ Amount by which total shipments of group z parts 

JZ3.U 

from plant a to consuming location(s) u during period j 
falls short of the minimum amount specified as desired in 
the sourcing rules 

[0059] ip Equations or Formulation Objective Function: 

[0060] Minimize: 

IIEX™W%~ ^YLYLsuBc jnaa L jmm ^YLYL^T^ + 

j m a e j m n a j m a v 

j "i a jzau jzau 

j m k q 



[° 061 ] Subject to: 

[0062] Sourcing Constraints: 



ma q 
ez 

l,(T Jmau +?,F jma J-S j2im <MAXPCT jzau H jzu 

m q 
? 

m q 

a. 

[0063] Capacity Constraints: 

YLCAPREQ jmaew P jmae <CAPACITY jaw 

m e 

[0064] Backorder Constraints: 

B jmkq = B (j _ X)mkq + DEMAND jnkq - £ F jmakq 

a 

[0065] Material Balance Constraints: 

Im = lu-»** + KECEIPT,mc + £ Z *7£ZZW*/W +£ Lj™ + 

" Z 2 QTYPERjmatnPlme 
.^component 

[0066] Non-Negativity Constraints: 

all AT, >y > 0 , where X is a generic decision variable and i, j etc. represent generic subscripts 

[0067] | n the supply chain linear programming model shown 

above, one problem that occurs is related to the need to 
improve runtime for large-scale production planning. Lin- 
ear programming models that represent division level 
production planning problems typically have variables and 



constraint equations measured in millions. Such large- 
scale production planning problems can pose serious 
problems in terms of runtime. Some supply chain plan- 
ning problems are naturally decomposable. For example, 
different manufactured parts may have independent bills— 
of-material. In these cases, supply chain models can be 
broken into multiple smaller independent linear program- 
ming models that can be solved much faster. However, 
many manufacturing industries assembly involve the use 
of common commodity components for multiple assem- 
blies (e.g., mechanical parts such as motors, cables, con- 
troller cards). These commodity items tend to link other- 
wise independent product groups and inhibit decompos- 
ing the resulting linear program. Furthermore, when sup- 
pliers provide multiple commodity items there may be 
contracts in place that specify contractual purchasing 
obligations for the aggregate of parts. Thus, otherwise in- 
dependent commodity items are linked by these contracts 
and similarly inhibit decomposition of the resulting linear 
program. 

[0068] There are existing methods for decomposing general (not 
necessarily supply-chain) linear programming models 
which are loosely linked (e.g. Dantzig Wolf decomposition, 



Benders composition "Linear Programming and Exten- 
sions", G. Dantzig). However, while these methods can 
improve runtimes they can be slow to converge and hence 
the improvements are often limited. The present method 
is an approach which uses a two stage linear program- 
ming method that specifically exploits the structure of 
supply-chain optimization models. 

[0069] | n many manufacturing industries (e.g., hard drives and 
servers in the computer industry) many components are 
used to build the final assembly, however, several of the 
components have unconstrained supply after some point 
in time (i.e., after the lead time required to obtain them). 
In a linear programming model of the supply-chain, these 
components act to prevent parallelization. The inventive 
decomposition method temporarily removes portions of 
the linear programming model (via relaxation of certain 
valid constraints) and solves a smaller relaxed version of 
the problem with improved runtime. By way of a two- 
stage algorithm (described below) certain constraints can 
be temporarily relaxed and then added back into the 
problem to obtain a solution to the complete model with a 
substantial improvement in runtime. 

[0070] The invention provides substantial runtime improvement. 



Large-scale production planning problems such as those 
at the central division planning level can pose serious 
problems in terms of runtime. The faster the application is 
able to run, the larger the number of scenarios that plan- 
ners are able to review in order to determine an optimal 
division build plan to commit. For some production mod- 
els instances, the invention has resulted in order of mag- 
nitude improvements in runtime (e.g., for one production 
model it reduced runtime from 13,201 to 933 seconds on 
an IBM Regatta p690). 
[0071] Figure 2 provides a graphical illustration of the inventive 
method. Figure 2 illustrates some simple bill-of-material 
(BOM) relationships in which assemblies, sub-assemblies 
etc. span four levels (L1-L4) of BOM and the overall linear 
program LP 22. Relaxation of material balance constraints 
and sourcing constraints for some of the L4 part numbers 
20 results in a linear program that is separable into three 
parts, LP1, LP2, LP3, which can be solved in parallel. Once 
these individual linear programs are solved, the basic 
variables in LP1, LP2, and LP3 can be fixed. "Basic vari- 
ables" is a term of art in the field of linear programming 
and describes those variables within the linear program 
model that are not set to zero as the linear program 



model is developed. Subsequently, the material balance 
constraints and sourcing constraints 20 relating to some 
BOM level L4 can be reinstated and the overall linear pro- 
gram can be solved (with most basic variables fixed) by 
leveraging the existing solution to the separable linear 
programs LP1, LP2, LP3. 

[0072] The following steps outline in detail the above-described 
method for decomposing a large-scale linear program- 
ming model for a supply chain optimization problem. 
More specifically, in item 300 the invention evaluates "ma- 
terial balance"and "sourcing"constraints associated with 
time dependent stocking points (PN/PLANT/PERIOD com- 
binations) to determine which ones can be temporarily re- 
laxed. The stocking point criteria relates to time depen- 
dent stocking points that include part numbers, locations 
of parts identified by the part numbers, and the time peri- 
ods when the parts will be available. 

[0073] The "material balance" constraints relate to changes in in- 
ventory of a specific part number at a specific location 
over time and the "sourcing" constraints maintain diversity 
in the source from which component parts are obtained. 
While these two constraints are used in this invention, the 
invention could use other similar constraints in different 



types of linear programs. For example, the constraints 
could involve any of the following characteristics dis- 
cussed below with respect to the criteria for choosing ap- 
propriate constraints. 

[0074] The criteria for choosing appropriate stocking points that 
are associated with the material balance and sourcing 
constraints includes selecting stocking points from the 
bottom of the bill-of-material. The bottom of the bill- 
of-material maintains the least complex components such 
as raw materials, unassembled components, etc., while 
the upper levels of the bill of materials include more com- 
plex part numbers such as partially completed sub- 
assemblies or complete assemblies. An additional crite- 
rion is whether the part number (PN) availability at the 
plant location (PLANT) is constrained by capacity con- 
straints (or similar constraints limiting availability). An- 
other criterion is whether the parts are available during 
the entire planning horizon. In other words, the PERIOD 
should be greater than (initial planning horizon date + 
shipping lead time + manufacturing cycle time) and the 
PN should be able to be built (i.e., is not inhibited) 
throughout the entire planning horizon. 

[0075] Thus, the invention identifies constraints that are associ- 



ated with the least complex parts within bills-of-materials 
used by the linear program as the constraints that can be 
temporarily relaxed. Additionally, this process identifies 
constraints that are associated with parts that have supply 
availability and lack capacity constraints as the constraints 
that can be temporarily relaxed. Further, this process 
identifies constraints that are associated with parts that 
are available during the planning horizon of the linear 
program as the constraints that can be temporarily re- 
laxed (the planning horizon includes an initial planning 
horizon, shipping lead time, and manufacturing cycle 
time). 

[0076] | n jtem 302 the invention relaxes the constraints that can 
be relaxed by setting the upper and lower bounds on the 
relaxable constraints to +infinity and infinity respectively. 
In item 304, the invention decomposes the resulting linear 
program 22, and solves the individual (smaller) linear pro- 
grams (LP1-LP3) in parallel on a multi-processor high 
performance computing server (e.g., simplex method). 

[0077] Next, in item 306, the invention replaces the basic vari- 
ables in the linear programs solved in item 302 with con- 
stants. More specifically, this processing selectively sets 
the lower and upper bounds for "basic" variables to their 



current value as solved in item 302. As mentioned above, 
"basic variables" is a term of art in the field of linear pro- 
gramming and describes those variables within the linear 
program model that are not set to zero as the linear pro- 
gram model is developed. In item 308, the invention puts 
the constraints relaxed in item 302 back into the linear 
program model and, in item 310, the invention solves the 
linear program (with the basic variables fixed in item 306) 
to determine the solution to the complete problem. 

[0078] with regard to the general flow of linear programming ap- 
plications illustrated in Figure 1, the current invention il- 
lustrated in Figure 3 relates to the solver illustrated in 
step 104. The method provides for improved computa- 
tional efficiency and opportunity for parallel processing of 
parts of the linear program. 

[0079] The present invention has been implemented on an IBM 
P690 server with AFS for data storage and using the AIX 
operating system. The steps for implementing the present 
invention are preferably programmed in C/C+ + . It should 
be understood by those of ordinary skill in the art, how- 
ever, that the representative invention is not limited to the 
above implementation and is independent of the com- 
puter/system architecture. Accordingly, the present in- 



vention may equally be implemented on other computing 
platforms, programming languages and operating sys- 
tems, and also may be hardwired into a circuit or other 
computational component. 
[0080] Thus, as shown above, the invention provides a method 
and system for solving a linear program having con- 
straints in a production planning system. The invention 
first determines which of the constraints can be tem- 
porarily relaxed based on stocking point criteria. The 
stocking point criteria relates to time dependent stocking 
points that include part numbers, locations of parts iden- 
tified by the part numbers, and the time periods when the 
parts will be available. The invention relaxes the con- 
straints that can be relaxed and decomposes the linear 
program into smaller independent linear programs. The 
invention initially solves the smaller independent linear 
programs with relaxed constraints (simultaneously in par- 
allel) to produce an initial solution. Next, the invention re- 
places variables in the linear program with constants 
based on this initial solution. After this, the invention re- 
stores the material balance and sourcing constraints and 
finally solves (re-solves) the linear program using the 
constants and with all constraints in place to obtain a 



complete solution of the linear program. 

[0081] The invention provides substantial runtime improvement. 
Large-scale production planning problems such as those 
at the central division planning level can pose serious 
problems in terms of runtime. The faster the application is 
able to run, the larger the number of scenarios that plan- 
ners are able to review in order to determine an optimal 
division build plan to commit. For some production mod- 
els instances, the invention has resulted in order of mag- 
nitude improvements in runtime (e.g., for one production 
model it reduced runtime from 13,201 to 933 seconds on 
an IBM Regatta p690). 

[0082] while the invention has been described in terms of the 
preferred embodiments, those skilled in the art will rec- 
ognize that the invention can be practiced with modifica- 
tion within the spirit and scope of the appended claims. 



